Equation of motion for dislocations with inertial effects 
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An approximate equation of motion is proposed for screw and edge dislocations, which accounts 
for retardation and for relativistic effects in the subsonic range. Good quantitative agreement is 
found, in accelerated or in decelerated regimes, with numerical results of a more fundamental nature. 

PACS numbers: 61.72.Bb, 61.72.Lk, 62.20.Fe 



I. INTRODUCTION 

Dislocation behavior in solids under dynamic condi- 
tions (e.g. shock loadingii^iS.) has recently attracted re- 
newed attention, ii^i^i^iiiSii^ partly due to new insights 
provided by molecular dynamics studies i^ii^i^ Whereas 
theoretical investigations mainly focused on the station- 
ary velocities that regular or twinning dislocations can 
attain as a function of the applied stress (possibly inter- 
sonic or even supersonic with respect to the longitudi- 
nal wave speed cl)^ i^°'^^ one other major concern is to 
establish an equation of motion^i^ '^^'^^ (EoM) suitable 
to instationary dislocation motions towards or from such 
high velocities, and which is computationally cheap. This 
would be an important step towards extending disloca- 
tion dynamics (DD) simulation o-'^^d^d^ to the domain of 
high strain rates, in order to better understand hardening 
processes in such conditions. 

The key to instationary motion of dislocations lies in 
the inertia arising from changes in their long-ranged dis- 
placement field, which accompany the motion. These re- 
tarded rearrangements take place at finite speed, through 
wave emission and propagation from the dislocation. As 
a consequence, dislocations possess an effective inertial 
mass,'^ which has bearings on the process of overcom- 
ing dynamically obstacles such as dipoles, etci^ d^d^ in- 
ertial effects are non-local in time, and are related to 
effective "viscous" losses. For small velocities where 
the EoM is linearj™ this relation takes the form of the 
Kramers-Kronig relations between the reactive and dissi- 
pative parts of the causal mass kerneh ^^d^d^ One major 
ingredient of the EoM should thus be the effective visco- 
inertial force exerted on the dislocation by its own wave 
emissioni^ii^ An EoM results from balancing it by the ap- 
plied stress, and by drags of various origins EoMs with 
effective masses, but which ignore retardation (e.g., Ref. 
3), cannot truly capture visco- inertial effects. Previous 
works on these questions having mainly been confined 
to the linear regime, their influence in the relativistic 
domain remains largely unexplored in spite of analytical 
progresses, partly due to the complexity of the formalism 
(especially for edge dislocations). 

Hereafter, Eshelby's EoM for screws with a rigid 
corCfi^ valid at small velocities, is first re-examined, and 



cast under a simple form which suggests a straightfor- 
ward regularization procedure for finite core effects. This 
allows us to appeal to previous results for point disloca- 
tions valid at high velocities.^ We then build in an heuris- 
tic way an EoM for accelerated or decelerated screw and 
edge dislocations in the drag-dominated subsonic regime, 
that consistently accounts for saturation effects at veloc- 
ities comparable to the sound speed. Results from the 
equation are compared to quasi-exact calculations from 
a numerical method of the phase- field type. Having in 
mind applications to DD simulations, the scope of the 
study is limited to continuum theory, so that dispersion 
effects due to the atomic lattice^ or to the influence of 
the Peierls potential,iS, are not explicitly considered. 



II. ESHELBY'S FORCE 

Within the Peierls-Nabarro model in isotropic 
elasticityf22,i^ and with the usual arctan ansatz for the 
relative displacement u{x, t) of the atoms on both sides of 
the glide plane, Eshelby computed the visco- inertial force 
F experienced by a screw dislocation of Burgers vector 
h, centered on position x{t) at time moving with a ve- 
locity V ~ i small compared to the shear wave speed 
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The dislocation is assumed to have a rigid core of half- 
width ^0- Then tg — 2Qo/cs is the time of shear wave 
propagation over the core width. The mass per unit dis- 
location length mo = /Lt5^/(47rcg) depends on the shear 
modulus ji. In Ref. [H (and in Ref. [3 as well), an incor- 
rect factor 1/2 is present in front of the second integral, 
and has been removed here. This factor is of no impor- 
tant physical consequence, save for different values of the 
linear response kernels; see below. 

That ([T]) is correct can be verified as follows. Starting 
from Eshelby's expression of the force as a double integral 
in Eq. (26) of Ref. [l^ . and expanding it to linear order in 
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the velocity w(r) or in x{t)—x{T), the following expression 
is easily obtained: 



Fit) = 2mo /' dr | ^=^| 



(At2+t2)5/2 



a;(t)-a;(r)]^ (2) 



where At = t — t. Using integrations by parts over r, 
each of Eq. Q and ^ can be put under the following 
irreducible form: 
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which shows them to coincide. 

By the same token, we check that ([T]) can be further 
simplified as: 
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By Fourier transforming Fit) [under the form ([3])] and 
by writing 

Fiuj) = [— CLi^m(ti;) — iujriiuj)] xiuj), 

we identify effective mass m(a;) and viscosity rjiuj) 
kernelsJ^ Their expression in closed form involves the 
modified Bessel and Struve functions Ki, Ii and C-i. 
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To leading orders in the pulsation to 
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where 7 is Euler's constant. Moreover, we observe that 

rii\uj\ ^ oo)/mo = 2/ts. (7) 

Result © coincides to leading order with Eshelby'Sfi^ as 
a; — > 0. The mass increase with wavelength as uj —> 
implies very different behaviors for, e.g., quasi-static and 
shock loading modes, since the latter involves a wider 
frequency range. We note that 77(0;) — > as w ^ 0, 
since losses should be absent from the model in the sta- 
tionary subsonic regimCf^^ The non-analytical behavior 
of the kernels at w = (due to \ui\), and its associated 
non-locality in time has been emphasized in Ref. 'T^. 

The finite "instantaneous" viscosity ([7]) stems from the 
first term in the R.H.S. of and is responsible for a 



velocity jump Av undergone by the dislocation when sub- 
jected to a jump AF in the applied force. ^'^^ From ([7]) 
we deduce: 



Av = 



AF tsAF _ ^^ CocsAF ^ 



riioo) 2mo 



The velocity jump ([5]) increases with core width. It was 
first predicted by Eshelby from his equation, ^nd can 
be understood as follows for a screw dislocation along the 
z axis: the force jump AF is due to a shear stress jump 
Aayz — AF/b attaining simultaneously all the points of 
the whole glide plane (e.g., as the result of shear loading 
applied on faces of the system containing the plane, paral- 
lel to the latter). Neglecting material inertia of the atoms 
on both sides of the dislocation plane, the medium under- 
goes an elastic strain jump Acjyz/ ^ = Asyz ^ Aw,„/cs, 
determined by a material velocity jump Avm- The latter 
is equilibrated through outward emission of a shear wave 
with velocity cg. On the other hand, the slope of the dis- 
placement function near the core is idu/dx) ~ 6/(2Co), 
so that Avm is related to the dislocation velocity jump 
Aw by Avm ~ Aw6/(2Co). Combining these relation- 
ships yields (jS]), up to a numerical constant factor. The 
same argument applies to other types of dislocations. In 
case of several relaxation waves (e.g., longitudinal and 
shear waves for an edge dislocation), that of lowest celer- 
ity controls the amplitude of the velocity jump. It should 
be borne in mind, however, that accounting for mate- 
rial inertia from the atoms on both sides of the glide 
plane results in an instantaneous inertial force of order 
Fi = 2mox to be added to By balancing the forces, 

it is seen that this force should spread the velocity jump 
over a short rise time 



At^ts. 



III. EQUATION OF MOTION 



(9) 



No expression analogous to ([T|) is available for edge dis- 
locations. However, Clifton and Markenscoff computed 
the force acting on a point screw or edge dislocation mov- 
ing with any subsonic velocity in an isotropic medium, 
that jumps instantaneously at instant t = t from rest 
to a constant velocity A generalization to anisotropic 
media is availablci^ To maintain its velocity constant, 
this dislocation must be subjected, at time t > t, to the 
time-decaying force 



it -T,V) = 



t - 7 



(10) 



where the function g depends on its character and on 
anisotropyi^ii^ We now construct heuristically a force 
for accelerated motion by interpreting such a motion as 
a succession of infinitesimal velocity jumps. Assuming 
that, for instationary motion, v in (jlOp can be inter- 
preted as w(t), the elementary force that would arise 
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from the elementary jump Sv(t) a,t t = t is: SF = 
[dF^M _ /dv{T)]6v{T) = g'{v{T))5v{T)/{t ~ r). 

Then, the total force experienced by the dislocation re- 
sults from integrating such elementary forces over past 
history: 



J — c 



dr 
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v{t). 
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Comparing pT|) to ([4]) shows, firstly, that the relevant 
"accelerations" at linear order are different. However, 
we remark that 2{d/dt){[x{t) - x{T)]/{t - r)} v{t) 
as i — > r, and moreover that for a screw dislocation, 
g'{v ~ 0) = moj^ Hence, since we interpret v in (jlOp as 
v{t), the numerator of the integrand in (jlip is correct 
at least for small velocities and for small times t ^ t. 
Its relevance for large velocities is demonstrated below 
through comparisons to full-field calculations. Next, in- 
tegral pTjl is singular at r = i, due the point-dislocation 
hypothesis at the root of (fTU|) . However, using (jH) as a 
physical motivation, we propose a regularization consist- 
ing in replacing the kernel 1/t in ifTTj) by + tg]^/^ 
where to, the counterpart of is in ©j is some time 
characteristic of sound propagation over a core diame- 
ter. In Sec. IIVI to is chosen alternatively proportional 
to ts = 2Co/cs and to ti, — 2(o/cl in the case of edge 
dislocations for illustrative purposes, whereas to is pro- 
portional to ts for screws. The proportionality factor, 
1/2 in all cases, is justified below. From a physical point 
of view, inertia is controlled by the slowest wave so that 
better results are expected using cs for all types of dis- 
locations. Given Eshelby's rigid-core hypothesis in 
and the approximations made, it would be pointless to 
refine this treatment. Another kind of regularization is 
used in Ref . [2^ (p. 195), which consists in replacing the 
upper bound t of integral ([TTjl by t ~ to (in Ref. [23, the 
integrand assumes that w ~ 0). 

With the above regularization the force eventually 
reads: 
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Its Fourier transform for small velocities where g'{v) ~ 
g'{0) yields, in terms of modified Bessel and Struve func- 
tions of order 0, 



m{uj)/g'{Q) = KoitoH) 
2e-'' 



(13a) 
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= ^|c.| + 0(tskP) 
77(|c^| ^ ^)/5'(0) = l/to. 
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The approximation therefore preserves the logarithmic 
character of the mass, and the viscosity, to leading or- 
der. The mass is slightly decreased, the constant mo/2 



in ([6a|) being absent. This difference is insignificant given 
the approximations made. In the limit of small velocity 
for a screw dislocation, our approximation amounts to re- 
taining in ^ the first integral only. In order to recover a 
correct velocity jump for screws, we must take to — ^s/2 
since the instantaneous viscosity p3cp is different from 
([7]). This "calibration" is used in the next section for 
screws and (somewhat arbitrarily) for edges as well. 

In the stationary limit, the visco-inertial force (|12p van- 
ishes. For V <^ Cs, the asymptotic velocity should be 
determined by a viscous drag force, mainly of phonon 
origin, ^'^ ^drago = Vov, where 770 is the viscosity. This 
force is modified (in the context of the Peierls-Nabarro 
model) by the relativistic contraction of the core, into 
-Fdrag(w) = ri{v)v. For subsonic velocities, ri{v) = 
rioD{0)/D{v), where;-iA 



D{v) ^ [A\v) + a'^iv/csy 



,1/2 



(14) 



with a = ?7oCo/(2'TioCs), is an effective viscosity- 
dependent core contraction factor, such that the 
core length in the laboratory frame reads: Q{v) — 
C,oD{v)/D{Q). The purely relativistic contraction factor 

A{v) is, with 7L,s = (1 - v^/clsY'^-^^ 

i lics/v)^ (47L - 7^^ - 27s - 7I) for edges, 
I ^7S lor screws. 

With this drag, and introducing the applied stress aa, 
the EoM finally reads: 



2tt 



dr- 



g'{v{T))v{T) 



+ Fd,,^ivit)) = baa, (15) 



[(t-r)2+i§]i/2 
where g{v) = 2mog{v), and wherei^ 

g{v) — {jg^ — l)/v, for screw dislocations, 

= (87L + 47L'-77s-67g-i+7s-3)ci/z;3 
-2[1 " ics/c1)'^]/v, for edge dislocations. 

This is our main result. By construction, it reproduces 
the asymptotic velocities of Ref. [ll|. 

We checked numerically that the replacement of Co by 
C(f ) in to does not change by more than a few percent the 
overall results described in the following section. Since 
this change in to would bring in nothing useful, we choose 
to use Co in ^0 in the following section. 



IV. APPLICATIONS 

Setting v{t) = Y,i Awj6'(i -ti), Eq. ^ is solved nu- 
merically for edge and screw dislocations, in an implicit 
way with a time step At = ti^i—ti small enough. Results 
are compared with numerical points obtained with the 
Peierls-Nabarro-Galerkin (PNG) approach^i^ used here 
as a benchmark. This method is less noisy than molec- 
ular dynamics, allows for full- field dynamic calculations 
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FIG. 1: Upper: Relationship between dimensionless applied 
stress a/fi and asymptotic velocity v/cs provided by the PNG 
code,— for an accelerated edge dislocation in the stationary 
regime of an accelerated dislocation (dots), compared to that 
predicted by Rosakis's Model I^^ (lines) for different viscosity 
parameters a, for a screw dislocation in the subsonic regime. 
Lower: Normalized velocity-dependent core width C('*')/C(0) 
measured under same conditions, cr is the Rayleigh velocity. 



of the displacement and stress fields in the whole sys- 
tem, accounting for wave propagation effects, and allows 
for better flexibility. We can thus, e.g., control bound- 
ary conditions by applying analytically computed forces, 
so as to prevent image dislocations from perturbing the 
simulation window. 

Firstly, to check the accuracy of the benchmark, 
asymptotic velocities of screw and edge dislocations were 
compared to the stationary predictions of Rosakis' Model 
iJi In the PNG method, the permanent lattice displace- 
ment field (which is part of the full atomic displacement, 
u) is relaxed by means of a Landau-Ginzburg equation, 
with viscosity parameter TypNG- An exact correspondence 
holds between this viscosity and Rosakis's viscosity pa- 
rameter a, namely rypNG = ctn/cs, as can be shown by 
specializing to one dimension the general field equations 
of Ref. 0. A 7-potential -f{u) = (1/2)70 sin^(7ru/6), with 
7o = (2/3)C446/7r^, is used. The material is an elasti- 
cally cubic material, with elastic moduli taken such that 
C44 = C12 = Cii/3 to insure isotropy. Due to the elas- 
tic correction made to the 7-surface potential in order 
to remove its quadratic elastic pait^^ the core at rest 
is a bit larger in the PNG results than in the Peierls- 
Nabarro solution. The time dependent core width C,(v{t)^ 
is measured from the numerical simulations by using 
6^/(27rC) = / dx[u'(a;)]^ (the value corresponding to a 




FIG. 2: Velocities vs. time for accelerated screw dislocations: 
white dots, PNG code; solid, Ref. 3 dots (in left curve only), 
Eq. lU; dash-dots, linear approximation to (I15|l : dashes, fully 
relativistic equation (|15p . 



core of the arctan type). Two-dimensional calculations 
are carried out using a simulation box of size 300 x 30 
6^, with a unique horizontal glide plane along Ox. Eight 
nodes per Burgers vector are used in both directions. 
Forces are applied on the top and bottom sides so as to 
induce shear on the unique glide plane. Free boundary 
conditions are used on sides normal to the Ox axis. Mea- 
surements are done near the center of the box, where the 
mirror attracting forces these sides generate on the dis- 
location, are negligible. The box is wide enough so that 
the dislocation accelerates and reaches its terminal ve- 
locity. Comparisons between PNG results and Rosakis' 
model are displayed in Fig. [T] for different viscosities a, 
in the case of an edge dislocation. The core scaling factor 
D{v) and the asymptotic velocity v/cg, are directly mea- 
sured from simulations under different applied stresses a, 
and compared to theoryjii The PNG asymptotic veloc- 
ities were found to be 5% systematically lower than the 
theoretical results. This correction is accounted for in 
the figure. The overall agreement is excellent. It is em- 
phasized that core contraction effects in the viscous drag 
[Eq. pi)) ] are required in order to obtain a good match. 

Next, comparisons in the accelerated regime are made 
with Eq. and with other models. Fig. [5] displays, 
as a function of time, the velocity of a screw disloca- 
tion accelerated from rest by a constant shear stress aa 
applied at t — 0. Low and high shear stresses are exam- 
ined. These stresses lead to terminal asymptotic veloc- 
ities Vt = v{t = 00) = 0.007 cs and 0.75 cs, computed 
from (fT5|) . The results displayed are obtained: (i) with 
the PNG approach (white dots); (ii) with Eq. using 
fully "relativistic" expressions of^(ti) and 13 (w) (dashes); 
(iii) with linear small- velocity approximations of g(w), 
but with the full expression of D{v), in order to empha- 
size the importance of relativistic effects in the retarded 
force (dash-dots, for the case vt = 0.75cs); (iv) with a 
previous EoMji using a typical cut-off radius R = 500 
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FIG. 3: Velocities vs. time for accelerated edge dislocations: 

white dots, PNG code; solid, Ref.Q; dash-dots, linear approxi- FIG. 4: Velocities vs. time for a decelerated screw disloca- 

mation to (|15p : dashes. fuUv relativistic equation (|15p . Curves tion. Comparison between the PNG method (white dots)^ 

obtained from (|15[l are duplicated, using either to computed and equation (|15p with (dashes) or without (dot-dash) fully 

with cs (upper), or with cl (lower), see text. relativistic expressions. 



nm in the logarithmic core term (solid) corresponding to 
a typical dislocation density of 10^^/m^. The result aris- 
ing from using ^ in the EoM is also displayed for the 
lowest speed (dots, left figure only). 

Figure [3] presents similar curves for an edge disloca- 
tion. For the latter, to is taken either as ts/2 or as 
— 2Co/cl, thus providing two limiting curves. 
The curves with tg provide the best matches, consistently 
with the above observation that the wave of lowest veloc- 
ity Cs < Cl should provide the main contribution to iner- 
tia. At low and high speeds, good agreement is obtained 
between PNG points and Eq. (fT5)) . provided that fully 
"relativistic" expressions are used for g{v) (especially for 
edge dislocations); otherwise, inertia is strongly underes- 
timated. In all the curves, the relativistic expression of 
the non-linear viscous terms was used. Moreover, varia- 
tions of the core width with velocity, implicitly present 
in PNG calculations, and ignored in the expression of 
to used in the visco-inertial term of (jlSp. are not cru- 
cial to accelerated or to decelerated motion (see Fig. U]) ; 
still, the core width shrinks by 20% during the accelera- 
tion towards Vt = 0.75 cl. On the other hand, retarda- 
tion effects in the effective mass are crucial: curves with 
non-local inertial forces are markedly different from the 
solid ones using the masses of Ref. |j, computed at con- 
stant velocity. The version of the PNG code used here 
does not include the above-mentioned effects of material 
inertia in the glide plane, so that the full-field velocity 
curves indeed display what resembles a velocity jump, 
like the EoM. Owing to ([5]), this lack of accuracy solely 
concerns the time interval between the time origin and 
the first data point: hence we can consider that the ve- 
locity jump is a genuine effect, and not an artefact, at 
least from the point of view of full-field calculations in 
continuum mechanics. However, we should add that, to 
our knowledge, this effect has not been reported so far 



in molecular dynamics simulations. Figure [3] displays the 
velocity of screw and edge dislocations decelerated from 
the initial velocity Vi. Comparisons between EoMs and 
PNG calculations are then harder to make than in the 
accelerated case. Indeed, the non-relativistic (resp. rel- 
ativistic) theoretical curves from Eq. (fT5| (dashed-dot) 
[resp. (dashed)] are obtained by assuming that an applied 
stress abruptly vanishes at t = 0. This induces a nega- 
tive velocity jump in the curves. This jump is larger if 
non-relativistic expressions are used, which demonstrates 
in passing the higher inertia (i.e. "mass") provided by 
relativistic expressions. The same loading was tried in 
the PNG calculations as well, but led to non-exploitable 
results due to multiple wave-propagation and refiection 
phenomena. Therefore, PNG curves for decelerated mo- 
tion were obtained instead using a somewhat artificial 
loading: the medium was split in a zone of constant 
stress, separated from a zone of zero stress by an immo- 
bile and sharp boundary. The dislocation is then made 
to accelerate in the zone of constant stress. Due to the 
finite core width, the boundary is crossed in a finite time 
~ C,/vi, which explains the smoothed decay of the veloc- 
ity in the PNG data points. This type of loading cannot 
be realistically implemented within the framework of Eq. 
(fT5)) because the dislocation core is not spatially resolved. 
Hence, though the curves strongly suggest that relativis- 
tic effects are as important in deceleration as in accel- 
eration, and that reproduces well the PNG points, 
the comparison between the latter and theoretical curves 
should be taken here with a grain of salt. On the other 
hand, the EoM of Ref. 4 (solid) is once again clearly im- 
precise. As a final remark, we expect our neglecting of 
retardation effects in the nonlinear viscous term of (fT5|) 
to induce an underestimation of damping effects. This 
may explain why the PNG curves decay faster than that 
from Eq. (fT5)) . 
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V. CONCLUDING REMARKS 

An empirical relativistic equation of motion for screw 
and edge dislocations, accounting for retardation effects 
in inertia, Eq. (fT5|) . has been proposed. We compared 
it, together with another available approximate EoM, 
to a quasi-exact numerical solution of a dynamical ex- 
tension of the Peierls-Nabarro model, provided by the 
Peierls-Nabarro Galerkin codei^ The latter was before- 
hand shown to reproduce very well the asymptotic veloc- 
ities of Rosakis's model 1^^ in the subsonic regime. The 
best matches with full-field results were found with our 
EoM, both for accelerated and for decelerated motion, 
thus illustrating quantitatively the importance of retar- 
dation and of relativistic effects in the dynamic motion 
of dislocations. To these effects, our EoM provides for 
the first time a satisfactory approximation for high ve- 
locities in the subsonic range. Our comparisons rule out 
the use of masses computed at constant velocity. One 



of the restrictions put forward by Eshelby to his EoM 
was its limitation to weakly accelerated motion, mainly 
due to the rigid core assumption.— Ours makes no at- 
tempt to explicitly overcome this simplification. How- 
ever, comparisons with full-field calculations, where the 
core structure is not imposed from the start, but emerges 
as the result of solving the evolution equation for the dis- 
placement field, shows that this rigid-core assumption is 
acceptable on a quantitative basis as far as inertia is con- 
cerned, at least for velocities high, but not too close to 
cs- 
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